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Neutron-rich gamma-ray burst flows: dynamics and particle 
creation in neutron - proton collisions 
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Abstract. We consider gamma-ray burst outflows with a substantial neutron component that are either dominated by thermal 
energy (fireballs) or by magnetic energy. In the latter case, we focus on the recently introduced AC' model which relies 
on magnetic reconnection to accelerate the flow and power the prompt emission. For both the fireball and the AC model, 
we investigate the dynamical importance of neutrons on the outflow. We study particle creation in inelastic neutron - proton 
collisions and find that in both models the resulting neutrino emission is too weak to be detectable. The inelastic collisions also 
produce y-rays, which create pairs in interactions with soft photons carried with the flow. In magnetically driven outflows, the 
energy of these pairs is radiated away as synchrotron emission. The bulk of the emission takes place at a few hundred keV, 
which makes it difficult to disentangle this signal from the prompt emission. In fireballs, however, pair cascading leads to the 
emission of y-rays with observer energy in the range of 2 - 20 GeV and a fluence well above the GLAST threshold. Therefore 
this emission can be a useful diagnostic of the nature of the outflow. 
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1. Introduction 

In recent years there has been significant progress in our under- 
standing of y-ray bursts (GRBs). The observational connection 
between supernovae and GRBs and studies of GRB host galax- 
ies provide compelling evidence for a connection between long 
GRBs and the death of massive stars (Van Paradijs et al. 2000; 
Woosley & Bloom 2006). The general scenario for long GRBs 
(for recent reviews, see Piran (2004); Meszaios (2006)) starts 
with core collapse of the massive star leading to the formation 
of a black hole surrounded by an accretion disk. The black hole 
- accretion disk system powers a developing outflow along the 
rotational axis, which accelerates to a bulk Lorentz factor of a 
few hundred, transferring its energy to the baryons contained 
in the flow. Dissipation of energy in the outflow leads to the 
prompt y-ray emission while the interaction of the outflow with 
the external medium results in the afterglow. 

The nature of the relativistic outflow is currently one of 
the most important open questions regarding GRBs. The high 
Lorentz factor, required to match the inferred energy density of 
the source and the observed non-thermal character of the emis- 
sion (the compactness problem; see e.g. Piran (2004)), implies 
that the ratio of energy to rest mass of the flow must be very 
high. In the widely used fireball model (Cavallo & Rees 1978; 
Goodman 1986; Paczyhski 1986) the outflow is a photon- 
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electron-positron plasma that is dominated by thermal energy 
and has a small baryonic load. Alternatively, the energy of the 
outflow may initially be dominated by Poynting flux (Usov 
1992). Such outflows occur naturally when a magnetized ac- 
cretion disk surrounds a black hole (Thompson 1994; Meszaros 
& Rees 1997; Spruit et al. 2001; Van Putten & Ostriker 2001; 
Vlahakis & Konigl 2001; Drenkhahn & Spruit 2002; Lyutikov 
& Blandford 2003; Lyutikov 2006; Uzdensky & MacFadyen 
2006). 

Neutrinos and y-rays may be useful probes to differentiate 
between fireballs and Poynting-flux dominated (PFD) outflows. 
The internal shocks that are believed to accelerate electrons in 
the fireball model will also accelerate protons to very high en- 
ergies, giving rise to neutrinos with energy >100 TeV through 
photopion production (Waxman & Bahcall 1997). In the ab- 
sence of a mechanism to accelerate protons to very high ener- 
gies these neutrinos are not expected in PFD outflows. In this 
paper we consider neutron - proton (np) collisions in neutron- 
rich flows and address the question whether neutrinos and y- 
rays created in these hadronic interactions can also be used to 
probe the nature of GRB outflows. 

GRB outflows are expected to be neutron-rich. In GRB 
central engines, the competition of positron capture on neu- 
trons and electron capture on protons favours a neutron-rich 
environment (Beloborodov 2003b; Pruet et al. 2003; Chen 
& Beloborodov 2007). Nucleosynthesizing interactions reduce 
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the number of free neutrons in the outflow, but a significant 
amount of neutrons remains in the flow until neutron decay 
becomes important (Beloborodov 2003b; Inoue et al. 2003). 
Deep in the outflow protons and neutrons are strongly coupled 
through nuclear scattering and behave as a single fluid that ac- 
celerates to high Lorentz factors. With increasing distance from 
the central engine the densities decrease until neutrons decou- 
ple and enter the coasting phase. Protons, being electromagnet- 
ically coupled to the flow, may be accelerated further. When the 
relative velocity between neutrons and protons is sufficiently 
high, inelastic np collisions are possible and lead to pion cre- 
ation. The pions decay into y-rays and neutrinos with observer 
energies in the ~10 - 100 GeV range. This mechanism has 
been investigated for fireballs (Derishev et al. 1999a; Bahcall & 
Meszaros 2000; Meszaros & Rees 2000; Belyanin et al. 2003; 
Razzaque & Meszaros 2006) but, to the best of our knowledge, 
not for PFD flows. 

The creation of secondary particles in inelastic np colli- 
sions can potentially be used to identify a substantial neu- 
tron component in GRB flows. Other ways to identify such 
a component that have been suggested in the literature are 
through signatures in the early afterglow of GRBs (Derishev 
et al. 1999b, Beloborodov 2003a; Fan, Zhang & Wei 2005), 
ultraviolet flashes generated in internal shocks in neutron-rich 
flows (Fan & Wei 2004), and observational signatures of a two- 
component jet that may be associated with neutron-rich MHD 
flows (Vlahakis et al. 2003; Peng et al. 2005). 

In this work we consider the 'AC' model as a specific model 
for PFD outflows. In this model the magnetic field configura- 
tion is similar to that produced by an inclined rotator (Coroniti 
1990; Lyubarsky & Kirk 2001) with field lines changing po- 
larity on a scale A ^ 2nc/Q., where Q denotes the angular fre- 
quency of the rotator. This model was recently discussed in 
connection to GRBs in a series of papers (Spruit et al. 2001; 
Drenkhahn 2002; Drenkhahn & Spruit 2002; Giannios & Spruit 
2005; Giannios 2006), where it was found that dissipation of 
the electromagnetic energy by magnetic reconnection can ac- 
count for both the bulk acceleration of the flow and for the 
prompt emission. 

The dynamics of fireballs and of outflows in the AC model 
are distinctively different. Fireballs are driven by radiation 
pressure and exhibit a period of rapid acceleration in which 
the Lorentz factor F oc r, where r denotes the distance from 
the central engine (Paczyhski 1986). The flow saturates either 
when there is no more energy available to further accelerate the 
baryons or when radiation and matter decouple at the Thomson 
photosphere. An analysis of the dynamics of neutron-rich fire- 
balls was recently presented by Rossi et al. (2006). The dynam- 
ics of neutrons in MHD flows was considered previously by 
Vlahakis et al. (2003) in the context of a different model for the 
outflow (Vlahakis & Konigl 2003) than the AC model consid- 
ered here. In the AC model, the acceleration of the flow is quite 
gradual and can be approximated with F oc r 1 ^ (Drenkhahn 
2002). Since acceleration of the flow is driven by magnetic 
forces, the flow can saturate far beyond the photosphere. It 
is expected that the difference in dynamics affects the number 
and the energy of secondary particles created in np collisions. 



Furthermore, the presence of a strong magnetic field can affect 
the interaction of secondary particles with the flow. 

Motivated by the fact that neutrinos and y-rays from inelas- 
tic np collisions could provide an indication about the nature of 
GRB outflows, we consider in this paper both fireballs and AC 
flows with a substantial neutron component. We investigate the 
dynamics of these flows and the creation of y-rays and neu- 
trinos in inelastic np collisions. In order to give an accurate 
comparison between the fireball model and the AC model, we 
consider both models here. Furthermore we use accurate fit- 
ting formulae for both the total and inelastic np cross sections, 
which has an important effect on the calculated fluences of sec- 
ondary particles. 

This paper is organized as follows. In section|2]we discuss 
the dynamical behavior of fireballs and of GRB outflows de- 
scribed by the AC model. In section[3]we consider particle cre- 
ation in inelastic np collisions. We discuss here the parameter 
space in which the mechanism is operational and we compute 
the fluences and energies of secondary neutrinos and y-rays. 
Detection prospects are discussed in section|4]and conclusions 
are presented in section [5] 

2. Dynamics of neutron-rich GRB flows 

Deep in the flow neutrons are strongly coupled to protons 
through elastic collisions, so that the two fluids behave as a sin- 
gle one. This np fluid is accelerated by conversion of thermal 
energy into kinetic energy in the fireball model and of magnetic 
energy into kinetic energy in the reconnection model. When the 
dynamical time of the flow becomes shorter than the np colli- 
sion time, the two fluids decouple and the neutrons enter the 
coasting phase. Provided that the flow has not already reached 
its terminal bulk Lorentz factor, the protons keep accelerating 
above the np decoupling radius, which results in relative mo- 
tion of the two fluids. 

The analysis of the effect of a neutron component on the 
dynamics is made separately for the fireball and the reconnec- 
tion model for the various stages of their evolution. Since the 
treatment of the mass flux is identical in both models, it is pre- 
sented first. 

2.1. Mass flux: protons and neutrons 

For an ultrarelativistic, steady, radial flow, assumed by both 
models under consideration, conservation of mass implies that 
the baryon outflow rate obeys 

M = M p +M B = 47i7- 2 mc(F p np+r n <) = 4nr 2 mcT p n' p {\ +£) , (1) 

where F p and F n stand for the bulk Lorentz factor of the protons 
and the neutrons, respectively, and n' p and n' n for their proper 
number densities. The masses of protons and neutrons are as- 
sumed equal m p ^ m n — m and £ stands for the neutron-to- 
proton mass flux ratio: 

M p F p n p 
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The ratio £ depends on the radius r since free neutrons decay 
into protons on a comoving timescale rg ~ 900 sec resulting in 

dM n = dM„ = _ 
dr r n cdf r n cT^ ' 

where t' stands for the comoving time. Taking into account 
that a proton is produced for every neutron that decays (i.e. 
dM n /dr = -dMp/dr), eqs. (f2| and 0} yield an expression for £ 
as a function of radius: 

dr T n CTj3 

From eq. (fl3 one can solve for the number density of protons 
and neutrons as a function of radius to find that 
1 M 



(4) 



1 + f 47ZT 2 mcT p 



and 



M 



(5) 



(6) 



11 1 + f Anr 2 mcY n 
The number density of the protons and neutrons is determined 
once their bulk Lorentz factor as a function of radius is derived. 
This is the topic of the next sections. 

2.2. The fireball 

In the fireball model most of the energy is initially stored in the 
form of thermal energy e, which is dominated by the energy 
density of radiation. The luminosity L of the flow is the sum of 
kinetic and radiation flux (e.g., Rossi et al. 2006): 



Anr 2 c 



r^(4e/3 + n' v mc 2 ) + Tln' n mc 2 



This expression can be rewritten as 
L = \nr 2 Y\cn' v mc 2 {\ + ^ + x) , 



(7) 



(8) 



where we have defined x = 4e/3n' p mc 2 . 

An important quantity for the evolution of the flow is the 
baryon loading parameter rj = L/Mc 2 » 1 where M (defined 
in eq. ([1])) includes both the contribution of the proton and the 
neutron fluid. Using expressions (Q} and ([8]l one derives the 
expression 

(i+£)77 = r p (i + x) + £r n . (9) 

Assuming that the flow starts from rest (i.e., r Pi o = r n> o =1) 
at an initial radius ro and initial neutron-to-proton ratio ft, the 
initial value for x is xq = 4eo/3n' pQ mc 2 — (1 + ft)(?7 - 1). 

As long as the flow is Thomson thick, radiation and par- 
ticles remain coupled and the evolution of the fireball is fully 
determined by the adiabatic lawQ 

e - e it) ' (10) 

P,0 



1 This expression does not take into account the increase of the pro- 
ton density due to neutron decay. The use of this expression is justified 
because, for the parameter space relevant for GRB flows, there is only 
a negligible fraction of neutrons that decays below the photosphere 
of fireballs. Hereafter, in the Thomson thick part of the flow, we set 
f = ft- 



From eqs. (0, (0 and ( TTOb one finds for the internal energy-to- 
proton rest mass ratio in the flow 



-*o(-^) 1/3 = (l+£>)(77-l)(^) 



1/3 



(ID 



Differentiating eq. (0 with respect to radius r and using eq. 
(fTTT i. one has an expression relating the bulk Lorentz factor of 
the proton and the neutron fluids in the optically thick part of 
the flow (see also Rossi et al. 2006) 



dT n 



2x 



3ft dr n 



(12) 



dr r 2x + 3 2x + 3 dr 

For the dynamics of the neutron-rich fireball to be fully de- 
termined, one needs to look closer at the momentum exchange 
between the neutron and the proton fluids because of np col- 
lisions. This has been studied by Derishev et al (1999a) and 
Rossi et al. (2006) who showed that when the two fluids have 
a relative velocity f3 K i, there is a drag force that accelerates the 
neutrons 



dr 



_E T 2 fl 2 

o 1 relPrel 



(13) 



where F re i ^ (F n /r p + F p /r n )/2 for ultrarelativistic flows and 
the total np scattering cross section a" to t is a function of jS re i. 
This expression accounts for the np interaction and does not 
depend on the acceleration mechanism (thermal or magnetic) 
of the flow. It can, thus, be applied to both fireballs and MHD 
flows. 

The np scattering cross section depends on the relative 
velocity of the two fluids. For np scatterings that take place 
with energies below the pion creation threshold, the scatter- 
ing cross section can with good accuracy be taken to scale as 
oc 1 / (ciySrei + C2y6j el ), while it remains almost constant for higher 
energies. The constants c\ and C2 are found by fitting to exper- 
imental data from Yao et al. (2006; see appendix A) 



cr tot = max 



0.19/^1 + 5.2^' 



(14) 



where & » 4 x 10~ 26 cm 2 . Our fitting formulae for cr tot are 
more accurate than the expressions used by Rossi et al. (2006), 
where the total np scattering cross section is substantially un- 
derestimated for y6 re iF re i < 1 (i.e. before np decoupling). This 
results in some differences in the dynamics close to the decou- 
pling radius. We find that the two fluids decouple over a nar- 
rower radial range (i.e. sharper decoupling). Furthermore, the 
fitting formula ( fT4b results in terminal neutron Lorentz factors 
that are ~10% higher than those found when we use the Rossi 
et al. (2006) expressions for the np scattering cross section. 

With eqs. ( fTTI ). ( fT2l and ( [TBI one has the complete descrip- 
tion of the dynamics of the fireball in the Thomson thick part of 
the flow (i.e. below the photosphere). In the optically thin part 
radiation and matter decouple and expression ( fTOb is no longer 
applicable. 

Since radiation pressure is the driving mechanism of accel- 
eration in the fireball, one would expect no further acceleration 
of the flow to take place above the photosphere. On the other 
hand, although most of the photons do not scatter with elec- 
trons above the photosphere, the electrons (outnumbered by the 



4 



Hylke Koers & Dimi trios Giannios: Neutron-rich GRB flows 



photons by a factor ~ 10 5 ) are still repeatedly scattered result- 
ing in a residual acceleration of the flow in the optically thin 
region. This residual acceleration is given by the expression 
(Beloborodov 2002; Rossi et al. 2006, appropriately modified 
to include the neutron fluid): 



dr n 



„dr n 



dr dr 



cr T L r 



167rFpV 2 mc 3 



1- 



vr p (r p h)r/ 



r n d^ 



1 + £ dr 



(15) 



where <x T is the Thomson cross section and L r = 16;tt 
stands for the radiative luminosity of the flow. The first term in 
the right hand side of the last expression accounts for the resid- 
ual acceleration from radiation and the second for the effect of 
neutron decay on the dynamics. Using eqs. and ©, we have 
for the radiative luminosity of the flow: 



L r = L\ 



/ _ r p + gr n x 



(16) 

The expressions ©, (fT~3T >, ( fT5l ) and ( fTol l describe the dynamics 
of the flow in the Thomson thin regime. 

2.3. The reconnection model 

In the magnetic reconnection model the flow is considered 
starting from the Alfven point r& and is dominated by Poynting 
flux. The luminosity of the flow is the sum of the kinetic and 
Poynting flux: 



Anr 2 c 



r 2 p (4e/3 + n' v mc 2 ) + r 2 n n' n mc 2 



+ c(rB) 2 



(17) 



where B is the magnetic field strength in the central engine 
frame, which is dominated by its toroidal component. 

A detailed investigation of the properties of a neutron-free 
flow is presented in Drenkhahn (2002) under the assumption 
of a cold flow (i.e. a flow where the term 4e/3 is neglected 
with respect the other terms in eq. (fTTTi). A full numerical in- 
vestigation showed that the dynamical description under the 
cold flow assumption is rather accurate (Drenkhahn & Spruit 
2002). Hereafter, we assume that the flow is cold. One should 
keep in mind, however, that though of moderate dynamical sig- 
nificance, the internal energy of the flow - dominated by the 
energy density of radiation - plays a crucial role for its pho- 
tospheric emission (Giannios 2006; Giannios & Spruit 2007). 
Furthermore, the cold flow assumption can overestimate the ac- 
celeration of the flow in the Thomson thin region by up to 50% 
in the limit that the internally dissipated energy does not stay in 
the flow but is efficiently radiated away (Drenkhahn & Spruit 
2002). More realistically only a fraction of the dissipated en- 
ergy is radiated away and the error we make in the Thomson 
thin region is smaller. 

Setting e = and using eqs. (Q} and ( fTTT i we have 



L = 



r P + ^r n 



Mc 2 + c(rB) 2 = 



r P + ^r n 



Mc\\ +o-), 



(18) 



where cr = (1 + £)(rB) 2 /(T p + ^T n )Mc is the magnetization 
parameter of the flow and stands for the Poynting-to-kinetic 
flux ratio. Using the last expression, the baryon loading of the 
flow is 



n = 



Mc 2 



1+f 



In the reconnection model, the flow is considered starting from 
the Alfven radius with magnetization cr () . In the inner part of the 
flow the very frequent np collisions lead to r p n — T n o = V^o- 
In view of eq. ([T9l >. one has that 77 = yfcroil + ctq) - cr^ 2 . 
The initial magnetization <x is, thus, an alternative means in 
parameterizing the baryon loading of a PFD flow. 

The radial dependence of the magnetic field strength is 
given by the induction equation that is appropriately modified 
to take into account the magnetic field dissipation through re- 
connection (Drenkhahn & Spruit 2002): 

d(rB) _ rB 
dr CTd 
Here, 



(20) 



2nTl 



0-+ 1 



(21) 



is the dissipation timescale of the magnetic field (in the central 
engine frame), Q stands for the angular frequency of the rotator, 
and s parameterizes the magnetic reconnection speed v lec . As 
in most models of magnetic reconnection, v rec scales with the 
Alfven speed v A , i.e. v rec = sva (see, for example, Lyubarsky 
2005). A nominal value used for e is 0.1. 

By combining eqs. (0J, ( fTTb . ( fT8b and d20l > one can elimi- 
nate the magnetic field B and derive an equation for the bulk 
Lorentz factor of the protons and the neutrons: 

dr p dT n 2 (..,-3/2 r , r \ . £p-£n d£ 

The last expression, in combination with eqs. and ( fTST l, de- 
scribes the the dynamics of neutron-rich flows in the reconnec- 
tion model. 



2.4. Results 

Having derived a closed system of equations that describe the 
dynamics of neutron-rich flows, we proceed with the investiga- 
tion of the dependence of their properties on the parameters of 
the flow for both fireballs and strongly magnetized flows. 

2.4.1. The fireball 

By numerically solving eqs. ( fTTT ). ( fT2b and ( fT~3l > in the Thomson 
thick part of the flow and eqs. ©, ( fjll ), ( fTBT ) and ( fTBT ) above the 
photosphere, one can follow the various stages of the neutron- 
rich fireball (see also Rossi et al. 2006). 

In Figs. (Q]) and ©, the bulk Lorentz factors of the proton 
and the neutron fluids are plotted as function of radius for dif- 
ferent values of the parameters of the fireball model. The latter 
are the luminosity of the flow L, the baryon loading 77, the ini- 
tial neutron-to-proton ratio £0 and initial radius ro of the flow. 
All the models studied have rj > 100 relevant for GRB flows. 

These low-baryon flows pass through an initial phase of 
rapid acceleration. During this phase, the neutron and proton 
fluids are strongly coupled and move practically with the same 
bulk Lorentz factor. Setting T p =t F n in eq. (fT2l we have 



dr n 



(19) = 
dr 



2x 



r 2x + 3(l +£0) 



(23) 
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Fig. 1. Bulk Lorentz factor of the protons (thick lines) and neu- 
trons (thin lines) for different values of the initial neutron-to- 
proton ratio ft an d radius ro of the fireball. At small radii, both 
protons and neutrons are in the linear acceleration regime (gray 
line). After np decoupling the neutrons saturate while protons 
can be further accelerated by radiation pressure. At r ~ 10 15 
cm the neutrons decay into protons that interact and decelerate 
the preexisting protons. 
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Fig. 2. Bulk Lorentz factor of the protons (thick lines) and neu- 
trons (thin lines) for different values of the baryon loading r\ of 
the fireball. For low baryon loading (high 77), the protons are 
accelerated to much higher bulk Lorentz factors than the neu- 
trons. For high 77, the saturation of the protons takes place close 
to the Thomson photosphere while the photospheric emission 
is very powerful. 



In the limit of x » 3(1 + ft)/2, radiation pressure leads to the 
well known linear acceleration of the flow as function of radius 
(cf. Goodman 1986; Paczyhski 1986; Piran et al. 1993): 

r p « r n = - , (24) 

ro 

If no n p decoupling were to take place, the bulk Lorentz factor 
of the flow would saturate at — n at the saturation radius 
r s = T]r . 

Note that although at small radii the numerical results fol- 
low the linear scaling d24l > closely, there are deviations from 
this scaling appearing for T p > 100 for the models presented 
in Figs. 1 and 2. Since eq. ( 1241 is exact for a fireball with a 
negligible number of baryons, finite-77 flows have bulk Lorentz 
factors T(r) < r/Vrj. 

At larger radii the density of the flow drops and np scat- 
terings become less frequent. When the comoving dynamical 
timescale becomes shorter than the np scattering timescale, 
the two fluids decouple and the neutrons are not accelerated 
any more. The relative velocity increases rapidly at decoupling. 
One can define the decoupling condition as r re i/3 re i = 1 . Setting 
this condition in ( fT3] l and using also ( 124b one finds for the de- 
coupling radius 

r np = 2.6 x I0^rl%^(^y 13 cm , (25) 
and for the Lorentz factor at decoupling 

F np = 2.6 x l^V^^P ■ 

If the flow reaches its terminal Lorentz factor before np decou- 
pling has taken place, both the neutron and proton flows coast 
with the same speed. 



For a flow with a sufficiently high 77, i.e. 

,>, cr , 360^V a y 4 (i^)" /4 , (27) 

the protons keep being accelerated after np decoupling has 
taken place while the neutrons coast with F n ~ F np . The bulk 
Lorentz factor at np decoupling F np provides a good estimate 
of the saturation Lorentz factor of the neutrons r n s . To quan- 
tify this statement, we have compared the analytical estimate 
for r np with the numerical values of r n at a large radius (here 
taken at r — 10 17 cm) and found that the two quantities agree 
with each other within ~25% for the (rather large) parameter 
space T] a - < i] < 3000, 0.01 < L 52 < 10, < ft < 10 and 
1 < r 0>7 < 100. 

When condition ( f2Tb is satisfied, the protons are further ac- 
celerated by radiation pressure after np decoupling until either 
all internal energy has been used or the flow crosses the pho- 
tosphere, where the flow becomes transparent with respect to 
Thomson scattering so that radiation and matter decouple. 

An estimate of the maximum Lorentz factor of the protons 
is given by assuming a neutron-free flow after np decoupling 
with luminosity L that does not include the kinetic energy of 
neutrons (i.e. L — L - F np ftMc 2 /(l + ft)) and mass flux M - 
M/(l + ft). The baryon loading of the decoupled proton flow 
is 

f)= — =*7(l+ft)-ftr n p. (28) 
Mc 2 

The acceleration of the proton fluid will saturate at 

r PiS = mint?), ?) rad ] , (29) 

where ?7 ra d = (Lcrx/47rromc 3 ) I/4 gives the terminal Lorentz fac- 
tor of the protons when the acceleration of the flow is limited 
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by photospheric crossing (Beloborodov 2002). This estimate 
takes into account the residual acceleration in the optically thin 
region discussed in section [2721 

At still larger radii of the order of rp - FnpCr^ ~ 10 15 - 10 16 
cm, neutron decay has an appreciable effect on the dynamics of 
the flow. The neutrons decay into protons and interact with the 
faster moving proton flow, thereby slowing it down. Note that 
at distances 10 17 cm, practically all the neutrons have decayed. 
The terminal Lorentz factor of the protons there is r p o < rj. 
For flows with fj > f) m d, most of the energy is not used to accel- 
erate the baryons (resulting in F pco <K rj) but instead appears as 
photospheric emission of the flow. 

Further out, the flow enters the afterglow phase where it de- 
celerates because of interaction with the circumburst medium. 
This last phase is not considered in this study. 

2.4.2. The reconnection model 

We now present the various phases of the development of the 
flow in the context of the reconnection model. The neutron-free 
flow has been studied by Drenkhahn (2002) and Drenkhahn & 
Spruit (2002). Here we focus on the dynamical effect of the 
neutrons. In Figs.[3]and[4l the bulk Lorentz factors of the pro- 
ton and the neutron fluids are plotted as function of radius for 
different values of the parameters of the reconnection model. 
These parameters are the luminosity of the flow L, the ini- 
tial magnetization <x () of the flow (that also parameterizes the 
baryon loading since 77 cr^ 2 ), the initial neutron-to-proton 
ratio £0 and the combination sQ. that parameterizes the recon- 
nection speed. 

The flow passes through an initial phase of acceleration 
where the neutron and proton fluids are strongly coupled and 
move practically with the same bulk Lorentz factor. Setting 
T p =i F n in eq. (1221 we have 



°0,2 =1 ' L 52 =1 



dTp 
dr 



sQJl 



r 1 3/2 



TTCF 2 



K 2 -r P ). 



(30) 



In the limit of F p <& <t 



3/2 



(i.e. the flow is still dominated by 
Poynting flux), the last equation can be integrated analytically 
to find (Drenkhahn 2002): 



3/2 vl/3 
r \ , 3/2) ' 



(31) 



The reconnection model predicts a gradual acceleration of the 
flow F ~ r 1 / 3 in the regime yfo\) <K T p <K o^J 2 with the bulk 
Lorentz factor of the flow given by 

r p - T n = (3eQ^ /2 r/7rc) 1/3 . 



(32) 



This expression is valid as long as the neutrons have not decou- 
pled from the protons and the flow has not reached its terminal 
Lorentz factor Too = a\J 2 at the saturation radius 



'0 



3eO 



'0 • 



(33) 



At larger radii the density of the flow drops and nuclear 
scatterings become less frequent. When the comoving dynami- 
cal timescale becomes shorter than the np scattering timescale, 
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Fig. 3. Bulk Lorentz factors of the protons (thick lines) and neu- 
trons (thin lines) for different values of the initial neutron-to- 
proton ratios fo and reconnection speed parameterized by eO 
in the reconnection model. At np decoupling radius the acceler- 
ation rate of the protons is enhanced. This effect is particularly 
pronounced for ^0 1 flows. At r ~ 10 14 - 10 15 cm, the neu- 
trons decay causing deceleration of the protons. 

S =l, m)=l, L 52 =l 



1000 r- 




r (cm) 

Fig. 4. Bulk Lorentz factors of the protons (thick lines) and neu- 
trons (thin lines) for different baryon loadings parameterized 
by the magnetization parameter <x in the reconnection model. 
The bulk Lorentz factor of the neutrons at np decoupling is es- 
sentially independent of cr , in agreement with the analytical 
estimate > 



the two fluids decouple and the neutrons are not accelerated 
any more. Their relative velocity yS re i increases rapidly around 
decoupling. As for fireballs, one can define the decoupling con- 
dition as r re iyS re i = 1 . Setting this condition in eq. ( [T3l and using 
also eq. d32l one finds for the decoupling radius 



r np =4.1Xl0^ (£ O)^cr-f(i±^) 



-3/5 



The bulk Lorentz factor of the flow at the decoupling is 

a/5/l+for" 5 



r np = iio^ 5 ( £ o)^(i±H)" 



(34) 



(35) 
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If the flow reaches its terminal Lorentz factor at r s before np 
decoupling has taken place, both the neutron and proton flow 
coast with the same speed. For a flow with a sufficiently high 
a"o, such that 



O"0 > <T C 



= 23L 2 . / 2 15 (eQ) 2 / 



-2/15 



(36) 



3/2 



the protons are further accelerated after np decoupling has 
taken place while the neutrons coast with T n ~ F np . The bulk 
Lorentz factor at np decoupling r np provides a good estimate 
of the saturation Lorentz factor of the neutrons F n s . Comparing 
the analytical estimate for F np with the numerical values of F n 
at large radii (taken here at r — 10 17 cm), we have found that 
the two quantities agree with each other within ~10% for the 
parameter space cr cr < cr < 300, 0.01 < L 52 < 10, < £ < 10 
and 0.01 < (sQ) 3 < 10. 

The critical value cr cr corresponds to baryon loading rj CI 
100. For baryon loadings 77 > 100 relevant for GRB 
flows, np decoupling takes place before the saturation radius 
has been crossed. In this case a substantial amount of magnetic 
energy is dissipated at radii r > r np , which is used to accelerate 
the protons. 

At the np decoupling radius the flow becomes effectively 
less baryon loaded and the protons increase their Lorentz factor 
more rapidly than the F p ~ r 1/3 scaling. This enhanced accel- 
eration is particularly pronounced in neutron dominated flows 
(where £0 » 1 ; see Fig.[3j. A similar enhancement in the accel- 
eration has been found by Vlahakis et al. (2003) in the context 
of a different MHD model for GRBs. 

Note that soon after np decoupling has taken place the flow 
crosses the Thomson photosphere. The protons keep accelerat- 
ing after the photospheric crossing in the magnetized flow since 
the acceleration is magnetic and not driven by radiation pres- 
sure as in the fireball model. At larger radii, the protons can 

3 /2 

reach bulk Lorentz factors in excess of the limit o~ that char- 
acterizes a pure proton flow (shown with dotted line in Fig. 3). 

At larger radii the neutrons undergo beta decay. For high-^o 
flows, at radius r ~ r np crg/£o the number of neutrons that have 
decayed is comparable with the initial number of protons in the 
flow and the effect of neutron decay on the bulk motion of the 
protons becomes appreciable. 

After magnetic dissipation has ceased and most of the neu- 
trons have decayed, all the available energy has been trans- 
ferred to the protons. The bulk Lorentz factor of the protons at 
large radii saturates to the value r Pi00 = cr^ 2 . This takes place 
at r ~ 10 cm. At these radii the flow is expected to enter the 
afterglow phase which is not considered here. 



3. Particle creation in inelastic neutron - proton 
collisions 

In the previous section we demonstrated that for low enough 
baryon loading, the neutrons decouple before the acceleration 
of the flow is completed in both fireballs and PFD flows. This 
leads to neutrons and protons developing relative motions and 
to energetic np collisions. Here, we study the production of pi- 
ons through inelastic np collisions in the relativistic outflow 
and the subsequent decay of pions into y-rays and neutrinos. 



We present analytical estimates for the secondary particle flu- 
ences and energies, and compare these estimates with numeri- 
cal results based on the model discussed in section|2] 

For the analytical estimates, we approximate the proton and 
neutron Lorentz factors as follows: 



F p s j ; F n =i min [r p , F np ] 



(37) 



where r np is the Lorentz factor of the flow at decoupling, p is 
a model parameter that allows us to consider the fireball model 
and the reconnection model together (p = 1 for fireballs and 
p = 1/3 for the reconnection model), and ro is a suitable length 
scale. For the fireball model ro is the initial radius where the 
fireball in injected, which is a free parameter of the model. In 
the reconnection model ro = 7Tc/3sQ.cr^ 2 is a length scale de- 
fined by the specific combination of the parameters - it has no 
deeper physical meaning but merely serves in rewriting the ex- 
pression (f32b in a more compact form. 

Using the unifying notation (|37| > for the bulk Lorentz factor 
of the protons and the neutrons, we express the np decoupling 
radius and the Lorentz factor at decoupling as: 



"( 



crLr 2 / \ 2n+] 



87rpmc J (l + ^)t] 



d-L 



np 



\2,f+l 

;/ 



(38) 



(39) 



v 87iro/wwc 3 (l + £)t] 1 
which combines eqs. d25l ). d26l i. ( l34t and d35l l. 

3.1. The pion production radius 

For sufficiently low baryon loading in the flow, pion creation 
in inelastic np collisions is possible after np decoupling and 
the subsequent acceleration of the protons with respect to the 
neutrons. We define the pion creation radius r n as the minimum 
radius where the relative velocity between decoupled neutrons 
and protons is large enough to create pions through inelastic np 
collisions. 

The production of a secondary particle with mass p re- 
quires center-of-mass energy yfs > 2mc 2 + pc 2 . Assuming that 
F p (r) » 1 and r n (r) » 1 at radii r > r np , and taking the np 
collision angle equal to zero (tail-on collisions), we express the 
center-of-mass energy V* as 



^ = mc 2 ( X 1/2 +X- 1/2 ) , 
where we introduce the useful quantity 

r P (r) 



X(r) = 



r»(r) 



(40) 



(41) 



From eq. ( f40b , we find that pions (which are the lightest 
mesons) can only be created if^(r) > x„, whereto =2.13 cor- 
responds to neutral pion production and;^± = 2.16 to charged 
pion production. We will use the average value % n = 2.15 in 
this work. Using the approximate proton and neutron Lorentz 
factors expressed in eqs. (l3"7T i. we find that 



Or — Xtt r np ! 



(42) 
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where the decoupling radius r np is given in eq. ( 1381 ). The radius 
from which pions can be created is thus substantially larger 
than the decoupling radius. Since the density of the flow and 
hence the number of np scatterings decrease rather steeply with 
radius, it is important to discriminate between r np and r„ when 
considering particle creation in inelastic np scatterings. 

Pion creation by np interactions occurs only when the pion 
creation radius r n is reached before the flow saturates. For the 
fireball model, saturation of the bulk Lorentz factor of the flow 
occurs either when there is no more energy available to fur- 
ther accelerate the baryons or when the flow crosses the photo- 
sphere (cf. eq. d29l>). It can be shown that saturation occurs be- 
yond the pion creation radius only if the baryon loading of the 
flow is sufficiently small. We express this condition as 77 > r\ n , 
where r\ n is the critical value for inelastic np collisions to oc- 
cur in the flow. Approximating the proton and neutron Lorentz 
factors with eqs. ( 137] ). we estimate that 

^-5.1xl0 2 4fr-yV(^o), (43) 
where 

0*6) = 0.85 {xn + £o) 3/4 (l + ft)" 1 (44) 

is a slowly-varying function normalized so that iff(l) = 1. In 
deriving eq. d43l we take neutron decoupling into account by 
using the neutron-free luminosity L and mass flux M as defined 
in section l2".4.1l The numerical investigation of the fireball dy- 
namics (see section [2~.4.1b shows that the proton Lorentz factors 
are substantially below the F p oc r scaling solution around np 
decoupling. As a result, the proton - neutron relative velocity is 
smaller and the pion production radius is pushed outward with 
respect to the analytical estimate (142b . This effect makes it more 
difficult to create pions in the flow and requires rj to be higher 
than the estimate d43l . Using the numerical model discussed in 
section |2] we find that inelastic np collisions in fireballs occur 
generally when 77/77,,- > 2, where r] K is expressed in eq. d43l . 
For neutron-rich flows (£0 ^ 3), neutron decoupling results in 
a relatively pure flow so that the protons follow the scaling ap- 
proximation d43l more closely and inelastic np collisions oc- 
cur already when 77/77^ > 1.5. Nevertheless, these results place 
quite stringent conditions on the fireball model parameters so 
that only a small fraction of GRB fireballs is expected to exhibit 
inelastic np collisions between bulk protons and neutrons. 

For the AC model we find that, similar to the fireball case, 
inelastic np collisions only occur for a sufficiently low baryon 
loading. We express this as cr > <x -n , where we use eqs. ( 1371 ) 
to estimate that 

/ 2 \ 2/15 

o-o* - 38 x L% 15 (eQ)f 5 ly-^-J . (45) 

We find that eq. (j45j is consistent with the critical value for 
(To obtained from numerical results on the proton and neutron 
dynamics (using the numerical model described in section |2]l. 

This value of <x „. corresponds to a critical baryon loading for 

3 11 

inelastic np collisions 77^ = crJ ~ 230 which is much lower 
than the critical value required in fireballs. Inelastic np colli- 
sions thus take place for a larger range of the parameter space 
in the reconnection model with respect to the fireball. 



The strength of any neutrino and y-ray emission that is a 
result of the decay of the products (mainly pions) of these col- 
lisions depends critically on the optical depth to inelastic np 
scattering. The calculation of this optical depth is the topic of 
the next section. 

3.2. Optical depth for inelastic np collisions 

The optical depth dr for a neutron with velocity c/3 n to scat- 
ter inelastically with a population of protons with velocity c/3 p 
and proper density n' p within r . . . r + dr is given by (see, e.g., 
Landau and Lifshitz 1971) 

dr = owrpn p (— jdr« — \x - -jdr, (46) 

where we assume in the last approximation that F p » 1 and 
r n » 1 and that the collisions are tail-on. 

At low center-of-mass energies the elastic and inelastic np 
cross sections are energy dependent. We find that for Xn ^ 
X ^ 10 (which is the range of interest here) the elastic cross 
section is well described with cr^ix) = 0.75o"/ ln^-, where 
<x = 4 X 10~ 26 cm 2 . A comparison between this approximation 
and experimental data on the elastic cross section taken from 
Yao et al. (2006) is presented in appendix|A] In the following, 
we express the inelastic np cross section as 

/ 0.75 \ 

O-inelCr > Xn) = & I 1 - — I . (47) 

We note here that the energy dependence of the np inelastic 
cross section has an important effect on the optical depth. If 
one assumes a constant cross section cr inel = 3 x 10~ 26 cm 2 
(as is often done in the literature) the optical depths are larger 
by a factor ~4 for both the fireball model and the reconnection 
model. Hence, the more realistic cross section adopted in this 
work leads to substantially lower estimates for the number of 
created particles. 

We consider, in general, the situation that neutrons coast 
with a constant Lorentz factor F n while protons are accelerated 
up to infinity with a Lorentz factor F p oc r p . Keeping p as a free 
parameter, we integrate eq. d46l ) through the flow to find that 

r{p) = jTd^l - ^j(x- 1 -X~ 3 )X- 1/P , (48) 

where eq. fl38T l was used to eliminate all parameters but p. We 
thus find that the optical depth for inelastic np scattering is in- 
dependent of any model parameters but the dynamical power- 
law index p. This result is valid for all outflows with F p oc r p 
and T n = const, provided that r s » r„. When r s > r„, such as 
in the fireball model, it represents an upper limit. 

For fireballs (p = 1) we find from eq. d48| > that t fb < 0.2, 
which is an upper limit because the flow saturates close to the 
decoupling radius. The situation is complicated by the fact that 
saturation of the fireball can be due to energy requirements or 
due to crossing of the photosphere. The numerical results pre- 
sented in section |2AT] indie ate that both effects cause the flow 
to accelerate considerably more slowly than the scaling approx- 
imation F p oc r near the pion creation radius r n . This pushes the 
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pion creation radius outward and decreases the optical depth 
for inelastic np scattering. We compute the optical depth nu- 
merically by a straightforward numerical integration of eq. ( |46b 
using the values of F p (r) and F n (r) obtained with the numerical 
mode0 discussed in section [2] We find that for the parameter 
space 3.5/7,, < 1 < 5000, 0.01 < L 52 < 10, 0.3 < £ < 10, 
and 1 < ro 7 < 100 the optical depth is approximated to within 
-25% by ' 



r i B ^ 0.11 fl - ^ 
»7 



(49) 



In the (rather favorable for frequent inelastic scatterings) case 
where rj = 5000 and £ = 5, the optical depth is r FB 0.1. 
For lower values of the baryon-loading parameter (T]/r]„ < 3.5) 
the optical depth is smaller than the value given in eq. j49l . A 
representative value for a fireball with 77 a few times the critical 
value r] n is t fb 0.05. 

For the reconnection model the saturation radius r s is typi- 
cally much larger than r„. We can therefore estimate the optical 
depth r AC for an inelastic np interaction assuming that the pro- 
tons are accelerated to infinity. (In principle this overestimates 
the interaction probability, but the difference is very small 
because the interaction probability decreases rapidly with r.) 
Inserting p = 1/3 in eq. (gHJl, we find that r AC * 8 x 10 3 . Fhis 
value is consistent with numerical results for flows with £q ~ 1 . 
For reference values of the parameters L52 = ^0 = o"o,2 = 
(e£2)3 = 1, we find also numerically that t ac = 8 x 10~ 3 . For 
high values of £ (neutron-rich flows), the extra acceleration of 
the flow after neutron decoupling (discussed in section 12.4.21 ) 
increases the optical depth by a factor few. We find that for 
the parameter space 1.5 < O"o/o"o )7r < 10, 0.01 < L52 < 100, 
0.3 < <f < 10, and 0.01 < (eQ) 3 < 10 the optical depth is 
approximated to within ~25% by 



T^O.Ol^ 



1/2 



(50) 



In particular, the optical depth increases to t ac 0.03 for 
very neutron-rich flows (£0 — 10). For 1 < <x /<x 0jr < 1.5 
pion creation is marginally possible and the optical depth is 
smaller than the value obtained by eq. (f50~b . For very pure flows 
(o"o/crrvr Z 10) neutrons decouple very early (before power- 
law acceleration r oc r p is reached), which results in an optical 
depth smaller by a factor ~2 than the estimate given in eq. (T50b . 

The obtained optical depth for inelastic np collisions is the 
first step in calculating the fluences of secondary pions and 
their decay products. The calculation of the fluences and ener- 
gies of stable decay products requires a model for the average 
number and average energy of neutrinos and y-rays created by 
np interactions. In the following sections we consider in detail 
the production of pions and the subsequent decay into neutrinos 
and y-rays. 



2 In the numerical analysis, we use a more accurate but also more 
elaborate approximation (see appendix [A) for the cross section than 
the one given in eq. l |47| (, which results in lower optical depths. 
Because dynamical effects, discussed in the text, have a larger influ- 
ence on the optical depth we use expression d47t > for simplicity to de- 
rive an analytical estimate. 



3.3. Pion production 

For the collisions studied in this work, the typical incident en- 
ergy of the proton measured in the rest frame of the neutron 
is /?p ~ 1 GeV/c. In this regime experimental data on pion cre- 
ation in np collisions is scarce and there is no unambiguous the- 
oretical framework. The available data (in particular, Prokoshin 
and Tiapkin 1957, Kleinschmidt et al. 1980, Daum et al. 2002; 
see also electronic data files available at the PPDS website 
http : //wwwppds . ihep . su : 8(5(51/) show a rapid rise in the 
single-pion cross sections just above threshold, and indicate 
that two-pion exclusive production cross sections are compara- 
ble to one-pion exclusive cross sections for incident proton en- 
ergies (as observed in the neutron rest frame K') p' p ~ 2 GeV/c. 
Furthermore, the ratio of 7r° : depends on energy; the ratio 
7r~ : 7r + is 1 : 1 under the assumption of nuclear isospin sym- 
metry. 

Based on the available data for incident proton energies 
p' p ~ 1 GeV/c we model the energy distribution and average 
number of pions resulting from anp collisions as follows. We 
estimate the ratio of created pions as n° : n + : if~ — 2 : 1 : 1 . 
Hence the average number of neutrino^ and y-rays resulting 
from a single inelastic np collision is: 



JV V = 1.0; N v . 



1.0; N v 



0.5. 



(51) 



Experimental data indicate that, for incident proton energies 
p' p = 1.14 GeV/c, the distribution of kinetic energy T = 
E - m n c 2 for n + mesons peaks around 0.6 r max , where T max 
is the maximum kinetic energy that can be carried by the pion 
(Kleinschmidt et al. 1980). For n° mesons in the center-of- 
mass (COM) frame this ratio is almost unity below p' p < 1 .06 
GeV/c and decreases to T pe; ± - 0.5 r max at p' p = 1.29 GeV/c 
(Prokoshin and Tiapkin 1957). Here we assume a constant frac- 
tion of 0.6 for all pion species and take the average pion energy 
for a single np collision in the COM frame K" equal to the 
peak energy: 



(O = <peak = °- 6 <max + 0Am n C 2 , 

where the maximum pion energy is equal to 
„ s — 4m 2 c 4 + m^c 4 

^r.max — « i— ' 

2 \S 



(52) 



(53) 



and eq. d40b relates the center-of-mass energy yfs tox- In these 
equations, the parameter^ provides the only reference to where 
the collision has occurred in the developing flow. 

When the angular distribution of pions in the COM frame is 
known, one can derive the full particle distribution of the decay 
products and transform this to the observer frame in order to 
find the secondary energy as observed on earth. However, there 
is to our knowledge no accurate parameterization of the angu- 
lar distribution of secondary pions created in np collisions. In 
the absence of such a parameterization we estimate the average 
observed energy of neutrinos by boosting to the observer frame 



3 Here and in the following denotes both muon-neutrinos and 
-antineutrinos (and similar for electron-neutrinos). 
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from an intermediate frame in which the secondary particles 
are assumed to be isotropic^ 

Pions are created approximately isotropically in the COM 
frame of the np collision. When neither pions nor their decay 
products are affected by the flow, as is the case for neutrino 
production in the fireball model, the distribution of the daugh- 
ter particles can be taken to be isotropic in the COM frame. 
In the AC model however, the strong magnetic field deflects 
the charged pions significantly since the pion gyration period 
is much shorter than the pion decay time. We assume that in 
this case the pions will be distributed isotropically in the frame 
comoving with the proton fluid. (Any randomized component 
of the magnetic field will further contribute to isotropization 
in this frame). Furthermore, in both the fireball model and the 
AC model y-rays from neutral pion decay will interact with 
the soft photon field of the flow, resulting in the emission of 
lower-energy photons. In the following sections these issues 
are discussed and estimates are presented for typical neutrino 
and y-ray energies. 

The decay of a charged pion also yields one electron or 
positron with energy ~35 MeV. These contribute to the y-ray 
emission which is discussed in section[ 

3.4. Observed neutrino energy 

In the fireball model pions do not interact significantly with the 
flow so that neutrinos from charged pion decay can be taken 
to be distributed isotropically in the COM frame K" of the np 
collision. For a given value of x, the average energy of the pro- 
duced neutrinos in this frame is then 



(54) 



where e° 35 MeV denotes the average neutrino energy in the 
rest frame of the decaying pion. The average observed energy 
(e v ) is obtained by boosting to the observer frame with Lorentz 
boost factor Fcom = yT p r n (appendixlBlcontains a summary 
of frames and Lorentz factors used in this work): 



Fcom«')^ r np x U2 (4')e° v 



1 + z m n c 2 



l+Z 



(55) 



where z is the redshift of the source. In the last equality, we 
approximate the proton and neutron dynamics by eqs. ( [37l >. 

For fireballs the flow saturates close to the pion creation 
radius and the bulk of the collisions occur when x is equal to 
the saturation value x$- For flows with large 77 > 700, where 
saturation is reached due to crossing of the photosphere, the 
terminal Lorentz factor of the flow equals F p s = ?) ra d (see eq. 



4 We note that, in the literature, there are various choices regard- 
ing the frame (e.g. the neutron rest frame or the proton rest frame) in 
which the energy distribution of y-rays and neutrinos is computed be- 
fore applying the final boost to the observer frame. Any intermediate 
frame leads to the same results in the observer frame provided that the 
angular structure of the particle distributions is taken into account. If 
an isotropic distribution is assumed, the choice of intermediate frame 
is important and depends on the physics. 



d29l>). In this case, the critical valuers can be estimated using 
the scaling law expressed in eqs. (l37l i: 



*s-4.0xL: 



-1/12,-1/12 1/3 
73 



52 



' 0.7 



1+f 



1/3 



(56) 



For lower values of 77, the saturation value Xs is smaller and 
reduces to the threshold values^ = Xn = 2.15 when 77 = 77^. 
Adopting the value x - Xs - 4 we find from eqs. ( 137] ). d5"2l 
and OBI that the neutrino energy in the observer frame can be 
expressed as 



(e v ) 



1 +z 



(57) 



where a™ 4.5 accounts for the non-zero kinetic energy of 
pions when they are created and for the fact that the particle 
distribution is not isotropic in the neutron rest frame. Using the 
same parameter range as in section [3T21 we find from a numeri- 
cal analysis that ff FB should be slightly higher than this estimate 
and we will adopt ff FB ^ 6 in the following. 

In the AC model the situation is more complex because 
charged pions interact with the flow before decay and because 
pions are created at various radii in the flow. Since the pion 
gyration time is much shorter than both the synchrotron cool- 
ing time and their lifetime, pions will isotropize in the frame 
K' comoving with the proton fluid without significant energy 
loss. In this frame, the secondary pions are injected with en- 
ergy {e' n ) = r^<<'>, where F p ' = ^/(2mc 2 ) is the Lorentz 
factor of the incident proton as observed in the COM frame. 
The observed neutrino energy is then given by the following 
expression: 

( ^ = g£^ >fr'"^-) «g, (58) 

v ' 1 + z m n c L 2(1 + z) 

where we approximate the proton and neutron dynamics by 
eqs. d37l i in the last equality. Note that the interaction with the 
flow results in a substantial increase in the observed energy of 
the secondary particles. 

For flows described by the AC model np collisions occur 
at various radii with different collision energies and different 
values for the Lorentz boost factor T. Therefore, we should av- 
erage the observed energy given in eq. ( 1581 over the developing 
outflow. We express the probability for an interaction to occur 
while ^ is in the range ^ . . .x + dx as T(x)dx- Since the scaling 
approximations d37] > describe the flow around decoupling quite 
well in the AC model, we use equation (08]) to estimate that 



(59) 



Averaging eq. (|58l over this distribution we find that the ob- 
served neutrino energy can be expressed as in eq. d57] i with 
a AC - 20. This is in good agreement with numerical results in 
the same parameter range as in section [3~!2l 



3.5. Reprocessing of y-rays: pair cascades versus 
synchrotron cooling 

While the flow is optically thin with respect to the emitted neu- 
trinos resulting from charged pion decay, this is not the case 
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for the y-ray photons that are produced by neutral pion decay. 
In the proton rest frame, the y-rays are injected with average 
energy (for a given value of x) 



' m n c L 

where = 70 MeV. Integrating over the developing flow as 
in the previous section, we express (e^) Be® and estimate 
analytically that B 3 for both the fireball and the AC model. 
This is consistent with numerical results. Hence y-rays have a 
typical energy 10/3 ~ 200 MeV in the proton rest frame and 
are ejected at radii r>r n \ not far from the Thomson photosphere 
of the flow. 

At these radii both fireballs and reconnection flows carry 
a soft photon field with characteristic comoving energy in the 
~1 keV range (see Derishev et al. 1999a and Giannios 2006 
for the fireball and reconnection model, respectively). Because 
of this intense soft photon field the flow is very optically thick 
with respect to these ~200 MeV photons, which are scattered 
and create pairs. In addition to the pairs that come from neutral 
pion decay, one energetic electron (or positron) is injected in 
the flow for every charged pion decay. We have included this 
contribution in the calculations that follow. 

In the fireball model the dominant cooling mechanism of 
the electron-positron pair is inverse Compton scattering. The 
upscattered soft photons create more pairs resulting in pair cas- 
cades. In the reconnection model the flow is dominated by 
Poynting flux and the energy density of the magnetic field is 
much higher than the radiation energy density (see also eq. (10) 
in Giannios 2006). As a result, the first generation of produced 
pairs cool down mainly through synchrotron emission. We dis- 
cuss the outcome of the y-ray injection separately for the two 
models. 

3.5.1 . Pair cascades in fireball 

Photons in the fireball are upscattered by pairs and absorbed by 
soft photons during the pair cascade. In the case of a saturated 
cascade, where all upscattered photons are absorbed, about 
~10% of the energy of the y-rays can be converted into rest 
mass of the pairs (Svensson 1987). More realistically the cas- 
cade is expected to be unsaturated, converting a few times less 
energy into rest mass of pairs (Derishev et al. 1999a; Belyanin 
et al. 2003). 

The result of these pair cascades is twofold. Each injected 
y-ray photon is reprocessed to multiple softer photons and the 
flow is loaded with pairs that contribute to its opacity. Although 
the saturation point of the cascade depends on the shape of the 
soft photon spectrum, we roughly estimate that photons with 
energies e' ~ 3 MeV in the proton rest frame are able to es- 
cape (Belyanin et al. 2003). The overall emitted spectrum will 
be broad and most energy is emitted with observer energies in 
the range e ~ r p , s (e'...10e')/(l + z) ~ (2...20)/(l + z) GeV. 
The strength of this component and its detection prospects are 
discussed in the next section. 

To estimate the importance of pair loading in the flow, one 
should compare the number of produced pairs with the number 



of electrons (or, equivalently, protons) pre-existing in the flow. 
The flow has £o neutrons per proton out of which a fraction t 
scatters inelastically. This results in %qt inelastic scatterings per 
proton. Every scattering results on average in ~ 1 y-ray photon 
(see eq. (IBTl i) with a typical energy 70/3 MeV in the proton rest 
frame. A fraction / ~ 3% of this energy is used in rest mass of 
pairs which results in ~ 10/3 f pairs per y-ray. 

By setting £o = 1 and using the values of f3 and t relevant 
for the reference values of the parameters for a fireball (see 
section l3~2l and the beginning of this section), one finds that the 
pair cascades result in ~0.15 pairs per proton. For the neutron- 
dominated case where £o = 5 we find significantly more pair 
loading, viz. ~2 pairs per proton. Note that we find significantly 
less pair loading of the flow because of pion decay compared 
to previous works. The main source for this discrepancy comes 
from the fact that, as we have shown in section [3~2l the optical 
depth for inelastic np scattering is about one order of magni- 
tude smaller than the value ~1 that is typically assumed in these 
studies. 

Most of these pairs are produced at large radii - and low 
densities - so that they do not annihilate but stay in the flow. 
For £o ~ 1> the contribution to the opacity from pair loading is 
at most moderate. On the other hand, for £o » 1 the number 
of produced pairs exceeds that of the pre-existing electrons. A 
fraction of those are produced below the Thomson photosphere 
and its location is pushed to larger radii. This can have some 
backreaction on the dynamics of fireballs that can be acceler- 
ated to slightly higher bulk Lorentz factors than those calcu- 
lated in section 12.4.11 where this effect is neglected. 

3.5.2. Synchrotron cooling in the magnetized flow 

We now turn our attention to the reconnection model. The typ- 
ical energy of the electron-positron pair produced by scatter- 
ing of a y-ray (resulting from neutral pion decay) with a soft 
photon is ~120 MeV which corresponds to a random electron 
Lorentz factor y e =s 200 - 300. The produced pair finds itself in 

a strongly magnetized flow with comoving B' a ^L/cr 2 Tl ~ 

10 6 G for typical values of the parameters and for the radii 
where most of the pion creation takes place. 

Under these conditions, the synchrotron cooling timescale 
of the pair t' s ~ 10~ 6 sec is much shorter than the Compton 
cooling timescale. The lack of pair cascades leads to negligible 
pair loading of the flow. The peak of the synchrotron emission 
is located at e' s = e%B'y\/m e c ~0.2...2 keV in the proton rest 
frame. At the radii where most of the pion production takes 
place, the bulk Lorentz factor of the protons is T p ~ 400 - 500 
which results in observer synchrotron peak in the sub-MeV en- 
ergy range. Keeping the rest of the parameters fixed to their 
reference values, we find that the synchrotron emission peaks 
at observer energy e s 120 keV for £o = 1 and at e s 600 keV 
for £o — 5. The spectrum is characteristic of fast (synchrotron) 
cooling particles with an exponential cutoff above the peak and 
a low-energy spectral slope of -1/2. The strength of this com- 
ponent and its detection prospects are given in the next section. 
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4. Detection prospects 

Using the results obtained in the previous section on the num- 
ber and energy of secondary neutrinos and y-rays created in in- 
elastic np interactions, we discuss the detection prospects here. 



4.1. Neutrinos 

We express the neutrino fluence as observed on earth as 

N n N v P Dp 



4nDj 



(61) 



where N v is the average number of neutrinos created per in- 
elastic np interaction (we add the contribution of muon- and 
electron-(anti)neutrinos given in eq. (f5Tb here), !P np is the in- 
elastic np interaction probability, D p is the proper distance, and 



N n 



ft 



1 + ft Tjmc 2 



3.3 x 10 



52 



_2ft_ 



£53*73 



(62) 



denotes the number of neutrons contained in the outflow. In the 
last equation, E denotes the total isotropic equivalent energy of 
the burst. Since !P np « 1 we express !P np r, where t denotes 
the optical depth for inelastic np collisions. 

We consider the optimistic case of a nearby energetic burst 
at redshift z =0.1. Assuming a universe that consists of matter 
and a cosmological constant, the proper distance D p is given 
by the following expression: 



D„ 



dz' 



Ho Jo V^A,o + ^ m ,o(l +z') 3 ' 



(63) 



where Qa,o and Q m o denotes the current density parameters 
of the cosmological constant and matter, respectively, and Hq 
is the Hubble parameter. Using the currently favored values 
n A ,o = 0.76, Q m>0 = 0.24, and H = 73 km s" 1 Mpc" 1 (Yao 
et al. 2006) we find a proper distance D p = 1.2 x 10 27 cm. 
Inserting this in eq. d68l l we find the following neutrino particle 
fluences for the two models: 



<D FB a 10" 



\0.05/\ 1 + ft 



$* c ^2x 10~ 5 



V0.01/V 1 + 



■ft 



£ 53O" 



3/2 



(64) 



(65) 



As discussed in section I3T21 a typical value for the inelastic np 
optical depth in the fireball model is t fb = 0.05 and for the 
reconnection model t ac = 0.01. The dependence on the model 
parameters, as obtained from a numerical analysis, is expressed 
in eqs. (|49j and (l50l . 

From eq. d57l i. the average neutrino energy as observed on 
earth is equal to 



(e v ) = 



np c v 



1 +Z 



(66) 



where a is a numerical factor that accounts for the non-zero 
kinetic energy of pions when they are created and for the fact 
that the particle distribution is not isotropic in the neutron rest 
frame while we boost with F np to the observer frame. Based on 



the results found in section [3~4l we take ff FB = 6 and a AC = 20 
for the fireball model and the AC model, respectively. Using 
eqs. ( f26b and d35l > for the Lorentz factors at decoupling we find 
that <e FB ) * 50 GeV and (e AC ) * 70 GeV for reference values 
of the parameters and a burst at redshift z = 0.1. These values 
depend only mildly on the parameters through the Lorentz fac- 
tor at decoupling r np but the value of a may change by a factor 
~2 depending on the burst parameters. 

Following Bahcall and Meszaros (2000) we estimate the 
number of interactions R Y in a large-volume neutrino detec- 
tor due to the diffuse background as R v = fyvKbO-yNt, where 
Kb = 10 3 ??i ; 3 denotes the burst rate per year, <x v = 5 x 
10 _39 (e v /l GeV) cm 2 is the neutrino interaction cross section 
and N, = l0 39 N t: 39 is the number of target protons in the de- 
tector. For reference values of the parameters and an average 
redshift z = 1 we find that R™ * 0.3 year" 1 and R AC 0.07 
year 1 for the fireball model and the AC model, respectively. 
Note that, for comparison with the literature, this estimate re- 
lies on the rather optimistic reference value of 1000 bursts per 
year leading to neutrinos through inelastic np collisions. 

The predicted diffuse neutrino detection rate for the fireball 
model is a factor ~5 smaller than the results found by Bahcall 
and Meszaros (2000). This is primarily due to the more accu- 
rate cross sections used in this work and the distinction between 
np decoupling radius and pion creation radius. This distinction 
also implies that the condition for inelastic np collisions to oc- 
cur (as expressed in eq. ( l43l ) is more stringent than the con- 
dition presented by Bahcall and Meszaros (2000). Therefore, 
the fraction of GRBs for which np decoupling occurs is ex- 
pected to be much smaller and the reference value Kb = 10 3 is 
not very realistic. For the reconnection model, we find that the 
expected neutrino fluence is typically lower than those for the 
fireball model by a factor ~5. This results from the fact that the 
pion production radius is much larger than the np decoupling 
radius, which is a very robust feature of this model. The condi- 
tion for inelastic np collisions as expressed in eq. (05]), on the 
other hand, is fulfilled in a large range of the parameters of the 
model. It is therefore expected that np decoupling occurs in a 
large fraction of GRBs for the reconnection model. 

4.2. Gamma rays 

Secondary y-rays resulting from np collisions are reprocessed 
by the flow due to interactions with the soft photon field (see 
section [3T5b . This results in pair cascades for fireballs and in 
electron synchrotron emission for AC outflows. The total en- 
ergy (in the frame of the progenitor) that is injected in the flow 
in the form of y-rays is equal to 



r p r p w n r«')e? 

— Y L np JV n t y t 



(67) 



which defines the factor y. We find that y ^ 0.5 for both the 
fireball model (for 77 ~ few r\ n ) and the AC model (for <xo ~ 
few ctq „). For reference values of the parameters this implies 
that the fraction of the burst energy that is converted to y-rays 
is roughly 5 x 10~ 3 for fireballs and roughly 2 x 10~ 3 for the AC 
model. We assume that the bulk of the energy given in eq. d6Tb 
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Fig. 5. Minimum value for the total isotropic burst energy Ei so 
for which the reprocessed y-ray emission (in the fireball model) 
is above the GLAST threshold, as a function of redshift. In this 
figure we have taken £o = foj =773 = 1 , and we have taken the 
burst duration equal to 10 s. 



mechanism can inject y-rays in the flow in a different region of 
the GRB parameter space. 

In the AC model the energy is radiated as synchrotron emis- 
sion with energy of a few hundred keV in the observer frame 
(see section [3.5.21 ). The corresponding energy fluence ~10~ 5 
erg cm -2 is lower than the expected prompt emission for a 
burst at z = 0.1 with the reference values adopted here and 
for a typical prompt emission radiative efficiency >0.1. This 
makes it very hard to disentangle this y-ray signal from the 
prompt emission. Of course this conclusion holds as long as 
the energy of the reprocessed y-rays is much less than the en- 
ergy of the prompt emission and the radiative efficiency for the 
prompt emission is larger than the energy fraction ~2 x 10~ 3 
transferred to y-rays by np collisions in the AC model. On the 
other hand, even though this synchrotron component is in gen- 
eral weak, it may have a substantial contribution to the prompt 
X-ray emission since its flux increases with decreasing energy 
as f v ~ v~ l/2 (i.e., following the characteristic slope of fast- 
cooling synchrotron emission). 



leaves the source after reprocessing, albeit in photons of lower 
energies. 

In section !?. 5. II we estimated that the y-ray emission from 
pair cascades in the fireball model is in the range 2-20 GeV 
in the frame of the progenitor. From this we estimate the y-ray 
number fluence O y from a source at proper distance D p as 



where e£ asc ^10 GeV is the average y-ray energy emitted by 
the pair cascades. For an energetic burst at z = 0.1 the num- 
ber fluence is cD r 10 3 crrT 2 which can be detected with the 
upcoming GLAST satellite that has an effective area ~10 4 cm 2 
at these energies (Gehrels & Michelson 1999). In fact, we find 
that this emission is detectable for a fairly large range of pa- 
rameters. In figure [5] we indicate, as a function of redshift, the 
minimum total isotropic burst energy for which the gamma-ray 
emission by this mechanism is detectable with GLAST. In pro- 
ducing this figure we have chosen reference values for the rel- 
evant model parameters and assumed a burst duration of 10 s. 

The isotropic equivalent energy carried by the prompt emis- 
sion at ~MeV energies of a typical GRB is in the range 10 52 - 
10 54 erg. This is only a lower limit for the isotropic equivalent 
energy of the ultrarelativistic flow which may well be a fac- 
tor ~10 larger than the energy carried by the prompt emission, 
depending on the unknown efficiency of the mechanism that 
generates the prompt emission. Therefore the minimum energy 
shown in figure [5] is not very restrictive and we expect that 
this emission is detectable for a fairly large fraction of GRBs 
in which protons and neutrons decouple. This conclusion also 
holds for high redshifts where the volume for GRBs to occur is 
largest. 

Apart from the collisions between bulk protons and neu- 
trons considered in this work, pions can also be created by 
nuclear collisions as a result of internal shocks in the sub- 
photospheric region of the flow (Meszaros & Rees 2000). This 



5. Conclusion 

In this work we have found that y-ray emission resulting from 
inelastic collisions between differentially streaming neutrons 
and protons and reprocessed by the flow may be a useful diag- 
nostic of the nature of GRB outflows. Provided that the baryon 
loading of the flow is sufficiently small, a few per mille of the 
burst energy is reinjected in the flow through np collisions in 
both the fireball model and in the AC model, which was used in 
this work as a specific model for GRB flows that are powered 
by magnetic reconnection. In the fireball model, the injection 
of these y-rays in the outflow leads to pair cascades and subse- 
quently to the emission of y-rays with observer energy in the 
range of 2 - 20 GeV /(l + z). In figure [5] we show the mini- 
mum total isotropic burst energy, as a function of redshift, for 
which this emission can be detected by GLAST. The constraint 
on the energy is not very restrictive and hence this y-ray emis- 
sion should be detectable for a fairly large fraction of the GRBs 
in which np decoupling occurs. In the AC model, synchrotron 
energy loss prevents pair cascading and the energy is radiated 
away at much lower observer energies of a few hundred keV. 
This component is expected to be dominated by the prompt y- 
ray emission. 

The neutrino particle fluence from n ± decay created in in- 
elastic np collisions in the fireball model is found to be an or- 
der of magnitude smaller than previous estimates. This is due 
to the more accurate cross sections for elastic and inelastic np 
scattering used in this work and the distinction between np de- 
coupling radius and the pion production radius. The neutrino 
fluence in the AC model is smaller by another factor ~5 due to 
the very gradual acceleration of the flow, which is a very robust 
feature of the model. The energy of neutrinos from np interac- 
tions in GRB outflows as observed on earth is in the range 50- 
70 GeV for reference values of the parameters, which is some- 
what higher than previous estimates. We find that the observed 
neutrino energy in the AC model is higher than in the fireball 
model because the strong magnetic field causes the charged pi- 
ons to isotropize in the proton rest frame rather than in the col- 
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lision COM frame. Unfortunately the neutrino emission in both 
models is so low that it is very difficult to use its properties to 
constrain the physics of GRB outflows. 

In both the fireball model and the AC model we find that 
inelastic np collisions occur only if the baryon loading is suffi- 
ciently low (see section [3~Tl >. For the fireball model, this condi- 
tion is quite restrictive and we expect that inelastic np collisions 
are only possible for exceptional bursts. On the other hand, in- 
elastic collisions occur for a large range of the parameters in 
the reconnection model. 

The above results rely on a proper understanding of the dy- 
namics of the flow. We have discussed the effect of neutrons 
on the dynamics of the flow in section [2] (some numerical re- 
sults are presented in Figs. [T||4j». We present a numerical model 
which includes the acceleration of the protons due to energy 
conversion in the flow, coupling of neutrons to protons by nu- 
clear scattering (and the dynamical decoupling of neutrons and 
protons) and neutron decay. To a first approximation the dy- 
namics of protons and neutrons can be described by the analyt- 
ical model given in eqs. ( l37l i. This model provides a useful esti- 
mate for the np inelastic optical depth (section [3~2l ) and the en- 
ergies of neutrinos (section [3~4"l i and y-rays (section [3~5l l. These 
estimates are generally in good agreement with results obtained 
from the numerical model described in section [2] (some differ- 
ences are discussed in the main text). The analytical estimates 
can be extended in a straightforward manner to any flow with 
T p oc r p and F n = const. 

From an observational point of view, the most promising 
conclusion of this work is that y-ray emission resulting from np 
interactions may provide a signature of the nature of the flow 
(section [3~5l l. The difference in energy of the reprocessed y-ray 
emission between the fireball model and the AC model results 
essentially from the difference in the ratio of magnetic energy 
density to radiation energy density. Therefore the energy of this 
emission appears to be a robust probe for the physics of GRB 
outflows. In this work we have estimated the y-ray energy and 
fluence for reference values of the burst parameters. A more 
detailed analysis is necessary to study the spectral properties of 
the emission and compare it with other emission mechanisms 
over a broad range of parameters. 

It was pointed out recently that a substantial neutron com- 
ponent in GRB flows may affect the properties of GRB after- 
glows (Belobororov 2003a). This provides a way of constrain- 
ing the physics of GRB outflows from afterglow observations. 
The numerical model discussed in this work can be used to 
study this possibility in more detail. Another interesting ques- 
tion is whether inhomogeneities in the flow can cause signifi- 
cant particle production through np collision in the AC model 
(for fireballs, this was discussed by Meszaros & Rees (2000)). 
These issues are left for future work. 
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Fig.A.l. Experimental data and approximations of the total, 
elastic and inelastic np cross sections. The thin lines show the 
approximations used for the analytical model; the thick lines 
show those used in the numerical computation. 



the COM frame K" of the np collision is moving with Lorentz 
factor 



COM 



(B.l) 



In the COM frame, protons and neutrons are moving in oppo- 
site directions with Lorentz factors 



2mc 2 



ill 

2\r n 



1/2 



2 r„ 



1/2 



(B.2) 



where we take the proton and neutron masses equal to m. In the 
main text we use K' to denote the rest frame of either the proton 
or the neutron. If K' denotes the proton rest frame, F p — 1 by 
definition and 



r; = I(.i + ^) = 2(r;') 2 -i. 
n 2\r„ r v j K n> 



(B.3) 



Appendix A: Cross section approximations 

In this work we use the following approximations for the total 
and elastic np cross sections: 



cr tot = max 



0.19^1 + 5.2^,' 



(A.l) 



0.75a- 

o- a (x> X k)= -j , (A.2) 

where cr = 4 x 10~ 26 cm 2 and^- = r p /r n . At energies below 
the pion production threshold (x < Xn) the elastic cross section 
o"ei = ""tot- Hence the inelastic cross section above the pion 
production threshold can be approximated with: 

/ 0.75 \ 

OfaelUf > Xn) = 0" t0 , - (T el = 0" 1 - . (A.3) 



In* 

In these equations, /3 re i and^ are related to the incident proton 
momentum in the neutron rest frame p' p as follows: 

A-ei = , Pp ; (A.4) 



p 2 + m 2 c 2 




(A.5) 



The approximations given in eqs. ( IA.1I ). iA.2i and (IA.3b are 
shown in figure IA.1I together with experimental data (Yao et 
al. 2006) and the approximation that was used to describe the 
inelastic cross section in the numerical analysis. 



Appendix B: Frames and Lorentz factors 

The Lorentz factor of protons and neutrons in the observer 
frame K are denoted with F p and T n , respectively, and we as- 
sume that both r p » 1 and F n » 1. In the observer frame, 



